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ABSTRACT 



Context. HD 150136 is one of the nearest systems harbouring an 03 star. Although this system was for a long time considered as 
binary, more recent investigations have suggested the possible existence of a third component. 

Aims. We present a detailed analysis of HD 150136 to confirm the triple nature of this system. In addition, we investigate the physical 
properties of the individual components of this system. 

Methods. We analysed high-resolution, high signal-to-noise data collected through multi-epoch runs spread over ten years. We applied 
a disentangling program to refine the radial velocities and to obtain the individual spectra of each star. With the radial velocities, 
we computed the orbital solution of the inner system, and we describe the main properties of the orbit of the outer star such as the 
preliminary mass ratio, the eccentricity, and the orbital-period range. With the individual spectra, we determined the stellar parameters 
of each star by means of the CMFGEN atmosphere code. 

Results. We offer clear evidence that HD 150136 is a triple system composed of an 03V((f ))-3.5V((f+)), an 05.5-6V((f)), and an 
06.5-7V((f)) star. The three stars are between 0-3 Myr old. We derive dynamical masses of about 64, 40, and 35 Mg for the primary, 
the secondary and the third components by assuming an inclination of 49° (sin^^ i = 0.43). It currently corresponds to one of the most 
massive systems in our galaxy. The third star moves with a period in the range of 2950 to 5500 d on an outer orbit with an eccentricity 
of at least 0.3. However, because of the long orbital period, our dataset is not sufficient to constrain the orbital solution of the tertiary 
component with high accuracy. 

Conclusions. We confirm the presence of a tertiary star in the spectrum of HD 150136 and show that it is physically bound to the 
inner binary system. This discovery makes HD 150136 the first confirmed triple system with an 03 primary star. 

Key words. Stars: early-type - Stars: binaries: spectroscopic - Stars: fundamental parameters - Stars: individual: HD 150136 
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1. Introduction 

NGC 6193 is a young open cluster surrounded by an H ii region 
in the AraOBl association. The core of the cluster hosts two 
bright and massive stars separated by about 10": HD 150135 and 
HD 150136. 

Several investigations in different wavelength domains (op- 
tical, radio, and X-rays) allowed to improve the knowledge of 
HD 150136. In the optical domain, Niemela & Gamen (200l) 
reported this star as a binary system composed of an 03 pri- 
mary and an 06 secondary. Indeed, spectral absorption lines of 
N IV /15203 and N v AA46Q3 - 19 (wavelengths are hereafter ex- 
pressed in A), as well as the emission line Niv /14058, were ob- 
served in the composite spectra. Th ese lines are conside red to 
be typical features for 02-03 stars dWalborn et al.ll2002h . This 
is a rather rare category in which only a few of these stars are 
reported as binary systems. Therefore, even though these ob- 
jects are considered, with the WNLh stars, as possibly belong- 
ing to the most massive objects, their actual masses are scarcely 
known. 

Located at a distance of 1.32 + 0.12 kpc jHerbst & HavlenI 

[1977), HD 150136 would thus harbour the nearest 03 star and 
would constitute a target of choice for investigating the funda- 
mental parameters of such a star Initially, the orbital period of 
the system was estimated to about 2.7 d whilst the mass ratio 
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between the primary and the secondar y stars w as measured to be 
equal t o 1.8 (unpublished, see Garm anv et alj ri980). The anal- 
ysis of Nieme la & Game 3 (2005) allowed the authors to refine 
these values to 2.662 d for the orbital period and 1.48 for the 
mass ratio. Moreover, the possibility of a third component bound 
to the system was envisaged (though no clear feature was de- 
tected) to explain that the He i lines in the spectra did not follow 
the orbital motion of any of the binary components. 

In the radio doma in, HD 150136 was r eported as a non - 
thermal radio emitter (iBenagUa et all 120061: fee Becked l2007h . 
which suggests the presence of relativistic electrons accelerated 
in shocks issued from a wind-collision region. T his non-thermal 
source could be of about 2" (i.e., the resolution. iBenaglia et al.l 
i2006 l) and could be associated with a clos e visual component 
separated from HD 150136 by about 1.6" jMason e t al. 199^, 
i.e., about 0.01 pc if we assume a distance of about 1.3 kpc. So 
far, there has however been no indication that this star is physi- 
cally bound to the binary system. 

The Chandra X-ray observations of the young open clus- 
ter NGC 6193 revealed that HD 150136 is one of the X-ray 
brightest O-type sta rs known (logLx ( erg s"') - 33.39 and 
Lx/iboi - 10"^"*, ISkinner et alJ l2005h . Its emission proba- 
bly results from a radiative colliding wind interaction and ap- 
pears to be slightly variable on a timescale smaller than one day. 
According to Skinner et al. (2005), the discrepancy between the 
orbital period and the timescale of the variation in this X-ray 
emission is probably an occultation effect either linked to the 
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Stellar rotation or to a putative third star with a rather short pe- 
riod. 

All these investigations thus suggest that HD 150136 is a 
multiple system rather than a binary one, even though evidence 
is still lacking. In this context, an intense spectroscopic mon- 
itoring at high resolution in the optical domain, first initiated 
to study the binary itself, was performed to bring new light on 
this system. We organize the present paper as follows. Section|2] 
describes the observations and the data reduction technique. 
Section [3] reveals the existence of a third star in the data of 
HD 150136 and focuses on the disentangling process in order to 
compute the individual spectra of the components. A first deter- 
mination of the orbital solution of the known short-period system 
based on an SB2 configuration is presented in Sect. |4] We also 
demonstrate that these parameters can be optimized by assuming 
the external influence of a third component. The possible orbit of 
this third star is addressed in Sect.|5]and an SB3 partial solution 
is derived. This discovery makes HD 150136 one of the first sys- 
tems where it is demonstrated that the third compo nent is phys- 
ically bound to the inner system (see also Trl6-104. lRauw et alj 
[2001, and HD 167971, Blomme et al. 2007). However, it is the 
first time that the primary is an O 3 s tar. We thus in S ect. [6l use 
the CMFGEN atmosphere code ( Hillier & Milled 1998h to derive 
the fundamental properties of the stars to improve our knowl- 
edge on the O 3 star parameters. Moreover, the configuration of 
the short-period binary seems to suggest a wind interaction zone 
possibly detectable through the He ii /14686 and Ha lines. We 
devote Sect.Qto the analysis of the variability observed in these 
lines by applying a Doppler tomography technique. Finally, we 
discuss the implications of our results in Sect.[8]and we provide 
the conclusions in Sect.|9] 



2. Observations and data reduction 

A multi-epoch campaign was devoted to observing HD 150136. 
Between 1999 and 2009, we collected 79 high-resolution, high 
signal-to-noise spectra of this star with the echelle Fiber-fed 
Extended Range Optical Spectrograph (FEROS) successively 
mounted at the ESO-1.5m (for spectra taken in 1999 and 2001) 
and the ESO/MPG-2.2m (for observations between 2002 and 
2009) telescopes at La Silla (Chile). This instrument has a re- 
solving power of about 48000 and the detector was a 2k x 4k 
EEV CCD with a pixel size of 15yum x 15fj.m. Exposure times 
were between 2 and 10 min to yield a signal-to-noise ratio over 
200, measured in the continuum on the region [4800 - 4825] A. 
FEROS is providing 39 orders allowing the entire optical wave- 
length domain, going from 3800 to 9200 A, to be covered. The 
journal of the observations is listed in Table [1] (including the in- 
dividual signal-to-noise ratios). For data reduction process, we 
used an improved version of the FER OS pipeline working under 
the MIDAS environment dSana et al.l l2006). The data normaliza- 
tion was performed by fitting polynomials of degree 4 - 5 to 
carefully chosen continuum windows. We mainly worked on the 
individual orders, but the regions around the Si iv /1/14089, 4116, 
Hen /14686, and Ha emission lines were normalized on the 
merged spectrum. We also note that the flatfields of the 2009 
campaign presented some oscillations of unknown origin. These 
oscillations represent about 3% to 5% of the continuum flux. We 
applied a filter to smooth these flatfields and thus to limit the 
effects of the undesirable oscillations on the observed spectra. 

Equivalent widths (EWs) and the Doppler shifts were es- 
timated on each spectrum by fitting Gaussian profiles on the 
spectral lines. For the computation of the radial velocities, we 



adopted the rest wavelengths from IConti et aP (Il977h . for the 
lines with a rest wavelength shorter than 5000 A, and lUnderhilll 
(1994) . for the other ones. 

3. Tliird component and disentangling 

On the basis of the high-resolution spectra analysed in the 
present paper, we report the unambiguous detection of a third 
star (T) signature for the system HD 150136 (Fig.[T]i. The data 
indeed exhibit a clear variation in the third component compared 
to the rest wavelength (He i /15875.62) shown by the dashed line. 
Moreover, the last observation confirms the periodicity of the 
variation, suggesting that the third star is also a member of a 
binary system. Although several analyses have suspected the ex- 
istence of this component, it has never been unveiled until now. 




5B60 5870 5880 5890 5900 

Wavelength [A] 



Fig. 1. Normalized spectra of HD 150136 zoomed on the [5860- 
5900] A region and vertically shifted for clarity. The He i /15876 
line displays a clear triple signature. Primary (?) and secondary 
(S) components are represented by filled and open cuxles, re- 
spectively. 

To separate the contribution of each star in the observed 
spectra, we apply a modifie d version of the disentangling pro- 
gram of ' Gonzalez & Levatol (1200 6) by adapting it to triple sys- 
tems. This method consists of an iterative procedure that al- 
ternatively uses the spectrum of a first component (shifted by 
its radial velocity) and then the spectrum of a second com- 
ponent (also shifted by its radial velocity) to successively re- 
move them from the observed spectra. The obtained result is a 
mean spectrum of the remaining third star. This program also 
allows us to recompute the radial velocities of each star with a 
cross-correlation technique, even at phases for which the lines 
are heavily blended. The cross-correlation windows are built on 
a common basis that gathers the positions of the He ii /14200, 
He I ^4471, Hen /I4542, Hei ^4713, Hei ^4921, Oiii /15592, 
Civ /1/15801 - 12, and Hei /15876 spectral lines. In addition, we 
add the Nv AA46Q3 - 19 lines for the primary star The radial 
velocities obtained with this technique are listed in Table [T] The 
resulting spectra are disentangled in a [4000 - 5900] A wave- 
length domain and at Ha. Though the He ii /14686 and Ha lines 
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Notes. RVp, RVs, and RVt are the radial velocities measured by cross-correlation, expressed in the heliocentric reference frame (no correction for 
the systemic velocities has been applied). We stress that the error bars on these values are about 5 km s"'. 



present particular variations along the orbital phase, we include 
them in the process. However, the assumption of constant lines 
(on which relies the disentangling method) is no longer satisfied, 
thus making these line profiles unreliable. 

To obtain the final individual spectra, we have to correct the 
disentangled spectra by the brightness ratio. For that purpose, 
we use the disentangled and the observed spectra to measure the 
EWs and constrain the spectral classification of the three compo- 
nents of HD 150136. This determinati on is based on the quantita- 
tive classification c riteria presented bv lConti & Alschuled (Il97lh 
and IContil(ll973bh . which uses the EW ratios of the Hei AAAll 
and Hen /14542 lines to define the spectral type. This leads to 
an 03-03.5 primary star, an 05.5-06 secondary star and an 
06.5-07 third star The weak Niii /1/14634-41 emission and 
the strong Hen /14686 absorption lines of the secondary and 
the tertiary components suggest adding the ((f)) suffix to their 
spectr al types. In addition, the criterion given by Wal born et al.l 
(I2OOI confirms the classification of the primary. Indeed, its 
spectrum is reminiscent of that of HD 64568, i.e., close to an 



03 V((r))^ However, the N iv /14058 fine has an intermediate EW 
(~ 160 mA) between the values measured on the Niii /i4634 (~ 
140 mA) and Nm AA6A\ {~ 250 mA) fines, thereby indicating 
a spectral type between 03V((r)) and 03.5V((f+)). The ((f*)) 
reports a spectrum with the Niv /14058 emission line stronger 
than the Nm /1/14634-41 lines and a weak Hen /14686 ab- 
sorption line whilst the ((f^)) refers to a spectrum with medium 
N m /1/14634 - 41 emission lines, the weak He n /14686 line, and 
the Si IV /1/14089 - 41 16 lines in emission. To compute the differ- 
ent brightness ratios required to determine the relative luminosi- 
ties of the stars, we compare the EWs measured on the observed 
spectra with the theoretical values issued from synthetic models 
whose spectral classification is the same as for our disentangled 
spectra. We stress that these models are computed from solar sur- 
face abundances (except for the nitrogen abundance). To avoid a 
possible bias due to differences between fits and models, we also 
focus on tables of IConti & Alschuleii (Il97lh . IContil (Il973ah . and 
IConti & LeepI (11974 1)" for the secondary and the tertiary compo- 
nents (the spectral type of the primary is not available in these 
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Fig. 2. Best-fit CMFGEN models (red) of the disentangled optical spectrum (black) of the primary component of HD 150136. 



tables). This comparison between the "observed" EWs and the 
"canonical" EWs for isolated single stars is performed on the he- 
lium and oxygen lines, yielding brightness ratios of Ip/ls - 2.54, 
h/h = 1.51, and Ip/h = 3.87. An error bar at 1 - cr of about 
0.14 is estimated between the "observed/canonical" ratios mea- 
sured from the different spectral lines. The individual spectra, 
corrected from the brightness ratios, are shown as black lines in 
Figs.Eltog] 

4. Inner short-period binary system 

With the radial velocities refined by the disentangling program, 
we compute the orbital solution of the short-period system 
formed by the primary (P) and secondary (S) components. We 
first re-estimate the orbital period of the short-term binary by 
using the Heck-M anfroid-Mers ch technique (herea fter HMM, 
iHeck et al.lll985L as revised bv iGosset et al.ll200lh on the ba- 
sis of RVs - RVp values. We obtain an orbital period of about 
2.67456 + 0.00001. We then use the Liege Orbital Solution 
Package (LOSFB to determine the SB2 orbital solution of the 



inner binary. Figure |5] displays the radial-velocity curves, and 
Table |2] lists the orbital parameters determined for an SB2 con- 
figuration. The systern shows no significant deviation (test of 
ILucv & Sweeney! 1 97 Ih from circularity, so we thus adopt e - 0. 
As a consequence, Tq represents the phase of the conjunction 
with the primary in front. We see that the dispersion of some ob- 
servational points around the theoretical radial-velocity curves is 
rather large, which is probably due to the influence of the third 
star 

Consequently, we compute the averages, for the primary and 
the secondary stars over each epoch, of the differences between 
the observations and the fitted radial-velocity curves. We also 
add to these values the constant systemic velocities quoted in 
Table |2] Moreover, we compare these results to the systemic 
velocities of the primary and the secondary components com- 
puted on each individual observing campaign. Both methods re- 
port comparable global trends of the velocities for the inner sys- 
tem. The adopted values are reported in Table[3] By plotting the 
mean systemic velocities of the inner system (see Table O and 



' LOSP is developed and maintained by H. Sana and is available 
atlittpV/www.science. uva.nl/~hsana/losp.html The algorithm is based 



on the generalization of the SB 1 m ethod of I Wolfe et alJ d 19671) to the 
SB2 c ase along the lines described in lRauw et^T32ob(ll) and lSana et al.l 
j2006h . 
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Fig. 3. Same as for Fig. |2]but for the secondary component of HD 150136. 
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Fig. 5. Radial-velocity curves corresponding to the short-period 
inner SB2 system as a function of the phase. Filled circles (open) 
represent the primary (secondary) component (see Sect. |4] for 
details). 



the mean radial velocities of the third component computed on 
each epoch as a function of time (Fig.|6]l, we detect a clear trend 
toward a periodical motion which is reminiscent of an additional 



orbit. Therefore, we suggest that HD 150136 is at least an SB3 
system. 



5. The outer system 

5. 1 . The outer system as an SB2 

We use the mean systemic velocities of the inner system deter- 
mined from each epoch as a primary component, and we con- 
sider the third component as the secondary companion. These 
velocities do not allow us to determine the real orbital period of 
the wide outer system with accuracy. However, we can estimate 
this value between roughly 2950 and 5500 d, i.e, between 8 and 
15 years. Consequently, computing an orbital solution is use- 
less without further high-resolution monitoring of HD 150136. 
Nonetheless, the measured velocities can provide information 
about the orbit, as well as on the mass ratio between the inner 
system and the third component. To estimate the mass ratio, we 
calculate the average between the systemic velocities of the pri- 
mary and the secondary stars and we take these values as radial 
velocities of the inner system (P-i-S). We focus on the radial ve- 
locities of both outer system components by fitting the following 



5 



L. Mahy et al.: Evidence for a physically bound third component in HD 150136 




4000 



4100 



4200 



4300 



4400 




y. 
3 




4400 



4450 



4500 



4550 



4600 



4650 



4700 



o 




4750 4800 4850 4900 4950 5000 




6550 



6560 



6570 



6580 




5500 



Fig. 4. Same as for Fig. |2]but for the third component of HD 150136. 
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Table 2. Orbital solution for the short-period binary system computed with LOSP and an SB3 program as described in Sects. |4] 
and|5] 



SB 2 solution 
Primary Secondary 



SB 3 solution 
Primary Secondary 



P [day] 


2.67454 + 0.00001 




2.67455 ± 0.00001 




e 


0.0 (fixed) 




0.0 (fixed) 




To [HJD-2450 000] 


1327.164 


± 0.003 




1327.161 ± 0.006 






1.615 ±0.021 




1.586 ±0.019 




y [km s"'] 


-19.7 ± 1.5 


-17.8 ± 


1.8 


-26.1 ±2.5 -24.4 ± 


2.5 


^:[kms-'] 


208.9 ± 1.7 


337.2 ± 


2.7 


211.2 ± 1.5 335.0 ± 


1.5 


a sin i [R^} 


11.0 ± 0.1 


17.8 ± 


0.1 


11.2 ±0.1 17.7 ± 


0.1 


M sin^ / [Mq] 


27.9 ± 0.5 


17.3 ± 


0.3 


27.7 ± 0.4 17.5 ± 


0.3 


rms [km s"'] 


13.15 




7.95 





Notes. The given errors correspond to 1 - cr. 



relation by a least-square method: 

RV^{4>) = ci 7;yp+s(0) + C2 

where RVp+s{ip), RVj(<p) are the systemic radial velocities of 
the inner system and the radial velocities of the third compo- 
nent, respectively, whereas ci - —^^^ with Mp+s and Mj 



the masses of the inner system and the third star, respectively. 
The linear regression gives us values of about ci - -2.60 and 
Q - -79.58 km s"'. From C2, we deduce yj ^ -20 km s"'. 
To verify the impact of the period on the value of the mass ra- 
tio, we have also computed several orbital solutions by assum- 
ing different orbital periods sampled between 2950 and 5500 d. 
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Table 3. Systemic velocities of the primary and secondary components, as well as the mean value over each epoch (for more 
explanations see Sect.|4|l. 



HJD 


#obs 


Primary 


Error Primary 


Secondary 


Error Secondary 


Mean 


-2450000 




[kms-'l 


[km s-'] 


[km s-'] 


[kms-'] 


[km s-'] 


1327.9031 


1 


-32.6 




-36.1 




-34.4 


2037.8496 


2 


-38.2 


12.5 


-49.1 


27.0 


-43.7 


2382.7319 


6 


-22.2 


4.6 


-24.0 


5.5 


-23.1 


2783.6601 


3 


-27.5 


2.2 


-23.9 


2.5 


-25.7 


3133.4644 


21 


-20.6 


1.4 


-16.0 


1.7 


-18.3 


3511.4421 


11 


-16.9 


0.8 


-13.8 


1.2 


-15.4 


3798.7490 


13 


-8.5 


1.2 


-2.2 


1.5 


-5.4 


3863.3497 


7 


-2.7 


1.3 


-5.2 


1.7 


-3.9 


5025.8661 


15 


-32.2 


1.2 


-33.6 


1.8 


-32.9 



Notes. The quoted errors correspond to 1 - cr. Second column represents the number of observations per epoch and the last five columns give the 
systemic velocity of the primary, the error on this value, the systemic velocity of the secondary, the error on this value, and the average on the 
primary and secondary systemic velocities, respectively. 



20M1 3000 
HJD-2 450 000 



Fig. 6. Radial velocities of the outer system represented as a 
function of the time. Filled circles/filled stars correspond to the 
systemic velocities of the primary/secondary component, whilst 
the open circles report the mean radial velocities of the third 
component (more details are provided in Sect. HJl. 



All these solutions led to a mass ratio of about 2.94. However, 
among the velocities exhibited in Fig.|6] the point corresponding 
to HJD = 2452 037.8496 (i.e., to the observing campaign made 
in 2001) displays a clear deviation probably due to the shortage 
of observations during this campaign. If we remove this point, 
we obtain by the linear-regression method a mass ratio of about 
3.1 and by the other method a mass ratio of about 3.2. We thus 
conclude on a mass ratio of about 2.9 + 0.3 between the inner 
system and the third star. Similarly, from this sample of orbital 
periods, we derive an eccentricity higher than 0.3 for the outer 
system. 



5.2. HD 150136 as a triple system 

To better constrain the orbital parameters of the inner system, 
we take the influence of the third component into account by 
fitting all the components together For that purpose, we deter- 
mine with another method the orbital parameters of the inner 
system. This technique provides the approximate orbital solu- 
tions for the two SB2 systems in HD 150136, i.e., the inner sys- 
tem composed of the primary and the secondary stars and the 
outer one represented by the mean systemic velocities charac- 
terizing the centre of mass of the inner system and the radial 



velocities of the third component. Therefore, we perform on all 
the three sets of data (P, S, and T) a global least square fit with 
an assumed circular orbit for the inner binary. For both systems, 
we follow the same procedure. We first deduce the approached 
orbital solution using the keplerian periodogram introduced by 
IZechmeister & Kiirsterl (2009). This method consists in a fit of a 
simple sine function in the true anomaly domain. The approxi- 
mate orbital solution is then globally refined using a Levenberg- 
Marquardt minimization. We search the period in [0.6r,2.0r] 
for the outer orbit where T is the total time span of the observa- 
tions. We use a constant step of frequency for both systems equal 
to 0.01 /r. Figure [T] exhibits the orbital solution for the short- 
term binary system by considering the influence of the third star 
We report in Table |2] the refined orbital parameters (SB3 solu- 
tion). Among these values, we notably observe a decrease in the 
error bars. Following an F-test on the ratio of the two of the 
fits (SB2 and SB3 solutions, see Table |2]i, we conclude that this 
solution is clearly improved in comparison with the one com- 
puted in Sect. |4] (see Fig. |5]l. The improvement is sufficiently 
marked to remain significant even though the F-test is not fully 
valid due to the non-Unearity. 




Fig. 7. Radial-velocity curves corresponding to the short-period 
inner system as a function of the phase according to the SB3 fit. 
Filled circles (open) represent the primary (secondary) compo- 
nent. The velocities are expressed in the frame of the centre of 
mass of the P+S system (see Sect.|5]l. 
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6. Atmosphere analysis of the three stars 

The three components of HD 150136 are then modelled sepa- 
rately by using the CMFGEN atmosphere code (iHilUer & Milled 
[1998). Non-LTE models including winds and line-blanketing are 
computed. The computations are done in spherical geometry and 
in the co-moving frame. A more detai led description of this code 
is provided in Hi llier & Miller! (Il998l) . For the velocity structure, 
we use TLUSTY hydrostatic structures in the photosphere as- 
sociated with a yS- velocity law v = Vco(l - R/r'f for the wind 
part (Voo being the wind terminal velocity). The TLUSTY mod- 
els are ta ken from the OSTAR2002 grid of models computed 
bv iLanz & Hubenv (2003) and = 0.8, typical for O dwarfs, 
is adopted (Repo lust et alj 12004) . The models include the fol- 
lowing chemical species: H, He, C, N, O, Ne, Mg, Si , S, Ar , 
Fe, and Ni with the solar composition of iGrevesse et al ] (l2007l) . 
except for N where the abundance is kept free. This is also the 
case for C and O abundances, but we see below that they are not 
well constrained and so are fixed to the solar values. We also use 
the superlevel approach to limit the required computation mem- 
ory. When the model atmosphere has converged, the synthetic 
spectrum is generated through the formal solution of the radia- 
tive transfer equation. For that purpose, a microturbulent veloc- 
ity varying linearly (with velocity) from 10 km s"', in the photo- 
sphere, to 0.1 X Voo, in the outer boundary, is used. The emerging 
spectra are compared to the disentangled ones, corrected for the 
brightness ratios, in order to constrain the main parameters more 
accurately. 

We thus determine the effective temperature (Teff), the \ogg, 
the luminosity, the N content, the projected rotational velocity, 
and the macroturbulence of each component with the following 
diagnostics: 

- the effective temperature: T^s is computed on the basis 
of the ratio between the Hei AAAll and Hen /i4542 
lines. However, other lines of helium such as Hei AAQ26, 
Hen /14200, Hei .14388, Hei ^4713, Hei AA92\, and 
Hei /15876 are also considered as indicators. Generally, we 
estimate the error in the models on the determination of that 
parameter to be equal to 1000 K. 

- the gravity: we use the wings of the Balmer lines (H(5, 
Hy, and HyS) as main diagnostics to constrain the \ogg. 
We refrain from using the Ho- line because it forms, at 
least partially, in the stellar wind. For the tertiary star, the 
determination of the gravity is more uncertain because 
of the disentangling process (see Sect. [8]l. Therefore, the 
uncertainty on the log g is equal to 0. 1 dex for the primary 
and the secondary and to 0.25 dex for the third component. 

- the luminosity: the luminosity is constrained by the V-band 
magnitude by adopting a distance of 1.32 kpc and a colour 
excess E{B-V) of 0.44 as explained in Sect.|3] The accuracy 
on the luminosity depends on the uncertainty on the distance 
of HD 150136. 

- the surface abundances: the nitrogen abundance is deter- 
mined on the N iii triplet in the 4500 - 4520 A range but 
also on the Niv /14058 and Niii ^1/14634-41 Une profiles. 
The typical uncertainty on the CNO surface abundances 
is of about 50%. To estimate these uncertainties, we run 
models with several different values of N and compare 
the corresponding spectra to the observed line profiles. 
When a clear discrepancy is seen on all lines of the same 
element, we adopt the corresponding abundance as the max- 



imum/minimum value of the possible chemical composition. 
Carbon abundance is estimated to be solar for the three stars, 
whilst oxygen abundance is not determined because of the 
lack of diagnostic lines. 

- the rotational and macroturbulent velocity: the synthetic 
spectra are successively convolved by a rotational profile 
and then by a Gaussian profile representing the macrotur- 
bulence. First, the projected rotational velocity i s computed 
with the method of I Simon-Diaz & Herrerol (12007 ) based on 
the Fourier transform. Then, the macroturbulent velocity is 
determined so that the shape of the He i AAl 1 3 line is repro- 
duced best. When this line is not visible in the spectrum, we 
use the Hei /i5876 and/or the Civ doublet /1/15801 - 12 as 
second indicator 

From y = 5.5 + 0.12 (iNiemela & Gamen"2005'), a mean 
{B — V) - 0.16 + 0.01, measured from Reed's catalogue (Ree^ 
2005), and a (B - V)q = -0.28 jMartins & Plezi 12006 ). we 
compute an absolute V magnitude for the entire system equal 
to My - -6.47;^Q3g by assuming a distance of about 1.32 + 
0.12 kpc. From the different brightness ratios, we then obtain 
Mv^ = -5.91!;^ ^|^, Mv, = -'^■90^_oil = -4-44!;|:^^ for 

the primary, the secondary and the tertiary stars, respectively. 
The bolometric correc tion is computed on the basis of the ex- 
pression given by Mar tins et al.l (|2005|) from a typical effective 
temperature determined as a function of the spectral type of the 
three stars (see Table|4]l. These values coupled with the absolute 
V magnitude of the three stars yield the bolometric luminosi- 
ties of about log (Lp/Lo) = 5M!_l \l log(Ls/Lo) = 5.32+^;|2, 
and log (Lj/Lq) - 5.01;^y [Q for the three stars. These luminosi- 
ties are in good agreement with the ones of main-sequence stars 
having similar spectral classifications. 

The best-fit model spectra of the three components are dis- 
played in red in Figs. |2]to|4l whilst the derived parameters are 
given in Table |4] However, the optical spectra do not give us 
enough information to constrain the wind properties of these 
components. We accordingly limit our investigation on the stel- 
lar parameters. Indeed, the Ho- and He ii /14686 lines are likely 
to be affected (to some extent) by the colliding-wind interac- 
tion zone, at least for the primary and the secondary compo- 
nents. Since the disentangling generates mean spectra on the en- 
tire dataset, the resulting Ha and He ii /14686 line profiles could 
therefore be different from those of single stars with the same 
spectral classifications. We stress that, as we see in Fig. |2] both 
lines display P Cygni profiles for the primary, which means that 
the stellar wind of this star could be rather strong in compari- 
son to its spectral classification. Therefore, we also investigate 
the UV domain, but the small number and the poor resolution of 
lUE data in the archives do not allow us to separate the individual 
contributions of the three components. Given the large ionizing 
flux generated by an 03 star, we thus assume that the flux col- 
lected in UV is mainly from the primary, whilst the He n /14686 
and Hff lines in the optical range serve to constrain the wind pa- 
rameters of the secondary and the third components. Therefore, 
we recommend considering the wind parameter values (M and 
Voo) listed in Table |4] as preliminary ones useful for computing 
the model atmospheres. 

7. Doppler tomography of the inner system 

The 03V((f )) - 03.5V((f+)) + 05.5-6V((f)) inner system with 
its small separation between both components likely favours in- 
teractions between the stellar winds. If we look closer at the 
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Table 4. Stellar parameters of the three components. 



Pa ra m p tp 


Pnmarv 


SppnnHarv 


L LI til y 


Sp. type 


03V((f ))/ 
03.5V((f+)) 


05.5-6V((f)) 


06.5-7V((f)) 


Teff [K] 


46500 ± 1000 


40000 + 1000 


36000 ± 1000 




5.86+°;? 


c ^9+0.12 
10 


5 01+"" 


R[Re] 


13.14!3«2 


9.54!|:^^ 


8.24+|-?5 


logg [cgs] 


4.00 + 0.10 


4.00 ±0.10 


3.50 ± 0.25 


Mspec [Mo] 


63.9!^^-^ 


n+27.7 


7.8+{2/ 




70.6!^'i* 


36.2![o 


27.0!3-0 


M sin^ / [Mq] 


27.7 ± 0.4 


17.5 ±0.3 




N/H[xlO-''] 


8.0±4.0 


3.0±1.5 


3.0±2.0 


V sim [km s '] 


171±20 


136±20 


72±10 


Vraac [km S"'] 


65±10 


41±10 


14±5 


M [Mo yr"'] 


1.0 X 10^' 


7.0 X 10^'* 


3.0 X 10-' 


v„ [km s"'] 


3500 ± 100 


2500 ± 100 


1200 ± 100 



Notes. The given errors correspond to 1 
nitrogen is N/H = 0.6 x 10"*. 



cr. The solar abundance for 



observed spectra, we clearly see that HD 150136 displays weak 
emissions in the He ii /14686 and Ho- lines. To better represent the 
emissivity of the inner system and notably the possible interac- 
tions between the wind of the primary and that of the secondary 
star, we use a Doppler tomography technique on those line pro- 
files. This technique maps the formation region of these lines in 
velocity space. The absorption line profiles of the primary and 
the secondary as well as the presence of the third component in 
the spectra do not simplify the representation of the emissivity 
in the inner system (see left panel in Fig.[8]l. Therefore, a prelim- 
inary treatment to eliminate these contributions in the observed 
spectra must first be applied. To realize this "cleaning" process, 
the disentangled spectrum of the third star, shifted by its radial 
velocity, is subtracted from all the spectra in our dataset. We 
also take the He ii /14686 line profiles of the synthetic spectra of 
the primary and the secondary obtained with the CMFGEN code 
(Sect. |6ll. We thus subtract them by shifting them with their re- 
spective radial velocities to remove the absorption parts of the 
lines. The resulting restored line profiles only take the emission 
part of the He ii /14686 line into account (Fig. |8] right panel). 
However, such a process must be used with caution, especially 
for the Ha Une. This line is indeed particularly dependent on the 
mass-loss rate and on the terminal velocity of the star Therefore, 
this cleaning process requires fully reliable absorption line pro- 
files, which is not the case for our three stars. The reconstructed 
emissions of the He ii /14686 and Ha lines are nearly symmetric 
and follow the orbital motion of the primary component. To anal- 
yse their entire profiles, we apply a Doppler tomography tech- 
nique to these two lines. 

The Doppler tomography technique assumes that the radial 
velocity of any gas flow is stationary in the rotating frame of 
reference of the system. The velocity of a gas parcel, as seen by 
the observer, can be expressed by a so-called "S-wave" function, 

v{(p) = -v^ cos(27T(f)) + Vy sin(2;7r0) + v~, 

where (p represents the orbital phase, (v,, v,,) are the velocity co- 
ordinates of the gas flow and is the systemic velocity of the 
emission line. This relation assumes an x-axis situated between 
the stars, from the primary to the secondary, and a y-axis point- 
ing in the same direction as the o rbital motion of t he secondary 
star (see e.g., Rauw et al. 2002; Linder et al.ll2008l) . Usually, the 
method of the Fourier-filtered back-projection based on Radon's 



transform is used in astrophysics. However, this technique has 
the disadvantage of generating "spikes" in the tomogram if the 
observations do not densely sample the entire orbital cycle. To 
limit these artefacts, we thus apply another Doppler tomography 
program based on the resolution of the matricial system Af - g 
where A is the projection matrix, / is the vector of emissivity, 
and § is a vector concatenating the observed spectra. The ele- 
ments of the projection matrix are equal to 1 if the projection of 
the S-wave onto the (v^, v^,) plane is verified or to otherwise. 
The solutions of the Af - g system are computed from an itera- 
tive proces s called SIRT (S imultaneous Iterative Reconstruction 
Technique. iKak & Slanevl l 20 01). This technique starts from an 
initialisation of the solution vector and computes a correction at 
each iteration that improves the solution of the previous itera- 
tion. These corrections are given by 



y(«+l) _ yl 



in) , 



■A'(g-A/«') 



where ju is a constant so that < ;U < 2 and || ■ || is the matri- 
cial norm. This method converges to a least-square solution but 
requires more computational time than the Fourier-filtere d back- 
proje ction (for more information on that technique, see 

Eonb. 

The tomograms are computed from multi-epoch spectra. 
Consequently, we have to correct the radial velocities to bring 
the spectra to the same systemic velocity on the entire campaign. 
The He ii /14686 resulting Doppler map (Fig.|9ll exhibits a single 
region of emission concentrated near the velocity of the centre 
of mass of the primary star The maximum of emissivity is lo- 
cated at (vj, V,.) ~ (50, -137) km s"'. For the Ha line, we see a 
similar pattern as for the He ii /14686 line. However, the analy- 
sis of the S-wave of the highest emission peak reveals a feature 
at (v'v, Vy) ~ (-24,-155) km s"'. These resulting tomograms 
do not show any structures reminiscent of a wind-wind interac- 
tion or of a Roche lobe overflow. Although the emission peaks 
are not clearly centred on the centre of mass of the primary, we 
can assume given the size of the stars that the emissions of the 
He II /14686 and Ha lines appear to be formed in an unperturbed 
part of the primary wind. This result coupled with the P Cygni 
profiles obtained from the disentangling thus suggests that the 
wind of the primary star could be relatively strong for its spec- 
tral type. We notice that t hese D oppler maps are close to those 
computed bv lRauw et akl (|2002|) for HDE 228766. 



8. Discussion 

8. 1 . The inner system 

The projected rotational velocities and the radii quoted in Table|4] 
suggest that both components of the inner system of HD 150136 
are in synchronous rotation, i.e., the angular velocity is simi- 
lar for both stars. Moreover, from the minimum and evolution- 
ary masses reported in Tables |2] and H] we derive an inclina- 
tion of about 47.2°+2 5 for the primary and of about 5l.3°tl'l 
for the secondary. Consequently, by assuming that our stellar 
parameters are close to the reality, we estimate the inclination 
of the inner system to be about 49° + 5°. This value favours a 
non-eclipsing system and agrees with the previo us inclination 
of about 50° given by iNiemela & Gamed (l2005h . We then fo- 
cus on the possible synchronous co-rotation of the inner system. 
From the orbital period and the radii, co-rotational velocities of 
about 276+42 ^ ' '^^'^ ^^'^-30 ^ ' the primary and the 
secondary components, respectively, are found. Such velocities 



require an inclination of about 38°!|^ and 47°!^^ to agree with 
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Fig.S. Left: 
HD 150136 



The observed He ii /14686 line profiles in the HD 150136 spectrum. Right: The restored He n A4686 line profiles in the 
spectrum when corrected for the absorption lines. 



the measured V sin/. These inclinations are close to the value es- 
timated from the masses. Therefore, we may conclude that the 
inner system is probably in synchronous co-rotation. 

By adopting a mass ratio of 1.6 and an a sin / = 28.9 Rq, we 
estimate the minimal radii of the Roche lobes of the primary an d 
the secondary from the expression given by Eggl etorJ (Il983l) . 
We find (TJrl sin /)p = 12.1 R^ and (/Jrl sin /)s = 9.8 Rq. If this 
result confirms that the secondary does not fill its Roche lobe, 
the case of the primary depends on the inclination of the system. 
By assuming an inclination of 49° +5° and by taking the radius 
given in Table |4] this component would not fill its Roche lobe 
either. 

8.2. The tertiary component 

Pr evious investiga t ions such as that of lNiemela & Garnenl (l2005h 
or iBenagha et al.l (12006) . mentioned the existence of a third 
component in HD 150136 . Sometimes, the close visual compan- 
ion located at about 1.6" (" Mason et al.lll99 8) was quoted as this 
tertiary star. However, from the present analysis, we conclude 
that this close visual star cannot be this third component re- 
vealed by our analysis. First of all, the spectral classification of 
the close visual companion is estimated to be later than B3, while 
we clearly report here an 06.5-7V((f)) classification. Moreover, 
if we assume an orbital period in the range of 8 to 15 years for 
the outer orbit, the possible a sin / separation between the third 
component and the inner system is estimated between 3500 and 
9100 Rq i.e., between about 16 and 43 AU. With such a sepa- 
ration and assuming a distance of 1.3 kpc for HD 150136, we 
calculate a typical angular separation (a cos /) between 14 and 
38 mas, much smaller than 1.6" as estimated fo r the close vi- 
sual companion reported by Mason et alJ (Il998l) . On the basis 
of the spatial separation associated with a mass ratio between 
the inner system and the tertiary component of about 2.9, this 
third component should be detec table by interfero metry accord- 
ing to the scheme given by Sana & EvansI (l20Tll) . An observa- 
tion made with SAM/NaCo (Nasmyth Adaptive Optics System, 
Paranal, Chile) in March 2011 allowed detecting a possible ob- 
ject at the limit of being significant (Sana et al. in preparation). 
This object would be located below 25 mas. Of course, this infor- 



mation must be taken with great care and needs to be confirmed 
by future observing campaigns. 

Although HD 150136 is now established as a triple system, it 
could still have a higher multiplicity. Indeed, the possible close 
visual component located at about 1 .6" could be also bound to 
the system. In order to determine whether HD 150136 is a dy- 
namically stable hi erarchical triple s ystem, we use the stability 
criterion quoted bv lTokovininI (l2004l) : 

^"0^(1 ~ ^out)^ > 5-Pin 

where fin and f out are the orbital periods of the inner and the 
outer systems, respectively, and gout is the eccentricity of the 
outer orbit. Even though we do not know precisely the exact val- 
ues of the parameters of the external orbit, we can still report, in 
the worst case, that this criterion is satisfied for the configuration 
of HD 150136. 

8.3. Evolutionary status 

The stellar parameters determined with the CMFGEN atmo- 
sphere code allow us to give the positions of the three stars in 
the Hertzsprung-Russell (HR) diagram. As shown in Fig.[TOl all 
these components are located relatively close to the ZAMS. The 
isochrones reveal an age of about 0-2 Myr for the primary and 
the secondary stars, whilst the tertiary would be slightly older 
with an age of about 1-3 Myr. Even though this discrepancy is 
barely significant, such a difference in age was already observed 
in several studies. Indeed, Martins et al. (2011, and references 
therein) reported that a possible reason to see this difference in 
the inferred age could be attributed to a different rotational rate 
between the components. By taking an inclination of 49° into 
account for the primary and the secondary component and by 
assuming that the tertiary is in the same orbital plane as the in- 
ner binary, we measure rotational velocities of about 227, 180, 
and 95 km s ' for the primary, the secondary, and the tertiary 
components, respectively. To represent them together, we use the 
evolutionary tracks computed with an initial rotational velocity 
of 300 km s"' (Fig. [Toll. However, evolutionary tracks based on 
the individual rotational rates are also investigated, but this does 
not change the conclusions: the discrepancy does not seem to 
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Fig. 9. Top: Doppler map of the Hen ^4686 line of HD 150136 
computed after removing the absorption line profiles of the pri- 
mary and secondary stars, as well as of the third component. 
The crosses correspond to the radial velocity amplitudes of the 
primary and secondary stars. The shape of the Roche lobe in 
velocity space (thick dashed line) was calculated for a mass ra- 
tio (primary/secondary) of 1.6. The Doppler map was computed 
with V; - km s"'. On the colourscale, red (blue) indicates the 
maximum (minimum) emissivity. Bottom: Same as for the upper 
panel, but for Ha. 

come from a different rotational rate. It has to be stressed that 
the evolution of the rotational velocity, hence the efficiency of 
rotational mixing, is probably quite different in a close binary 
system where the stars corotate synchronously, on the one hand, 
and in a rather isolated third star, on the other. 

This discrepancy in age is strengthened by determining the 
stellar parameters (see Sect. |6]l. Even though the parameters of 
the primary and the secondary stars agree with those of stars on 
the main sequence, we observe for the third component that its 
luminosity corresponds to a main-sequence star, whilst its grav- 
ity is more reminiscent of an evolved star However, log g is de- 
termined on the basis of the wings of Balmer lines, and its esti- 
mate could be wrong if the wings of these lines are poorly repro- 
duced by the disentangling process. In the case of HD 150136, 



we see that, due to the broadness of the Balmer lines, the spectral 
separation between the different components is too small to sam- 
ple the full width of these lines. Therefore, the deduced value of 
the log^ could be affected by the disentangling process. 
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Fig. 10. HR diagram reporting the position of the three compo- 
nents of HD 150136. The evolutionary tracks computed for an 
initial ro tational velocity of about 3 00 km s and the isochrones 
are from lMevnet & Maeded (l2003h . 

Moreover, the parameters of the third star also show an evo- 
lutionary mass that is higher than the spectroscopic one, even 
though both values are lower than what the mass ratio between 
the inner system and the third star seems to indicate. We notice 
that the spectroscopic mass is derived from the radius and the 
log g of the star Consequently, a wrong estimate of log g affects 
the determination of the spectroscopic mass. If we consider the 
evolutionary mass of the inner system (primary and secondary), 
we find a total mass of about 100 Mq. By adopting the mass ratio 
we have derived between the inner system and the third star, we 
obtain a mass of about 35 Mq for the tertiary component. This 
value is close to the evolutionary mass, which would show that 
the spectroscopic \ogg would be wrong. This dynamical mass 
would also lead to a log g value of about 4.1, which is more rem- 
iniscent of a main- sequence star. 

Finally, if we assume that the disentangling does not af- 
fect the true value of the logg, we have to look for the real 
phenomena that could expla in such a value. From the differ- 
ent evolutionary schemes of 'Eg gleton & Kiseleval ([1996'), it is 
however unlikely that the third star has a giant luminosity class. 
Indeed, according to these authors, a massive triple system such 
as HD 150136, i.e., with Mp > Ms > Mj, should first undergo 
a Roche lobe overflow in the inner binary. Since the primary 
star is more massive than the other ones, it should be the first 
to evolve, thus implying the mass transfer in the inner system. 
Therefore, the tertiary component could not be evolved, whilst 
the primary and the secondary still remain on the main-sequence 
band. Furthermore, the absence of modified surface chemical 
abundances and the low V sin/ derived for the third compo- 
nent also do not support an alternative scenario where the third 
component would have been captured by the inner system after 
having been ejected by a supernova kick or by dynamical inter- 
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actions in the neighbourhood of HD 150136. All in all, it thus 
appears that the tertiary component and the two other stars be- 
long to the main sequence and that the most probable cause of a 
wrong estimate of log g is the disentangling process. 

8.4. Colliding wind and non-thermal radio emission 

The presence of an 03V((f ))-03.5V((f+)) star in the inner sys- 
tem seems to be mainly responsible for the variations observed 
in the Hen /14686 and Ha lines. However, the proximity of the 
05.5-6V((f)) secondary could imply a wind-wind interaction 
zone between the two stars, even though the disentangling and 
the Doppler t omography render th at unlikely. According to the 
formalism of IStevens et al.l (Il992h and the parameters given in 
Table H] we expect that the powerful wind of the primary star 
crashes onto the surface of the secondary star This scenario is 
in line with the finding of ISkinner et al. I (l200l . The collision of 
the primary wind with the secondary surface could generate an 
X-ray emission associated with the secondary inner surface. The 
short timescale X-ray variability could then result from occul- 
tation effects very much like those observed for CPD-41° 7742 
jSana et al ] |2005E . Though this scenario seems to agree with the 
results of the Doppler tomography where we have seen that the 
emission is only due to the primary star, the position of the con- 
tact zone remains unknown because the wind parameters of the 
secondary star are uncertain. Therefore, we use the mass-loss 
rate of about M - lO^^Mg yr ' derived bv iHowarth & Prinial 
([1996') as a function of stellar luminosity, and we assume an 
inclination of 49° for the system. These assumptions lead to 
a probable radiative primary wind, with a cooling parameter 
X = 0.87. On the basis of differ ent values for the fundamental 
parameters, ISkinner et al. I (l200l derived = 0.2, also in agree- 
ment with a non-adiabatic regime for the X-ray emission from 
the inner wind-wind interaction region. On the other hand, any 
X-ray contribution coming from the colliding wind region lo- 
cated between the inner system and the third star should a priori 
be produced in the adiabatic regime, at least in a very large frac- 
tion of the long eccentric orbit. 

The confirmation that HD 150136 includes a third star is 
worth discussing in the context of its non-thermal radio emis- 
sion. The stellar wind material is quite opaque to radio photons, 
which are unlikely to escape if the absorbing column along the 
line-of-sight is too thick. The radius of the t - 1 radio photo- 
sphere can be est i mated on the basis of the formula given by 
IWright & Barlow! (Il975h an d Leitherer et al. ( 1995), following 
the same approach as i De Becker et al.j (i2004) in the case of the 
candidate binary system HD 1681 12. Using the stellar wind pa- 
rameters given in Table |4l we calculate the radio photosphere 
radii, respectively for the primary, the secondary, and the tertiary 
components. These radii are plotted as a function of wavelength 
in Fig. [TT] When discussing the radio emission from massive 
stars, one should keep in mind that stellar winds are radio emit- 
ters, but in the thermal regime. On the other hand, the wind-wind 
interaction in massive binaries could pro duce non-thermal radi- 
ation in the form of synchrotron radiation. iBenaglia et alJ (l2006t) 
reported on a clearly non-thermal spectral index between 3.6 
and 6 cm, but the spectral index at longer wavelengths (13 and 
20 cm) is more typical of a thermal spectrum. We therefore focus 
on the values derived at shorter wavelengths. While we are deal- 
ing with a triple system, let us separately consider two possible 
wind interaction regions: a first one within the short period sub- 
system (hereafter, SPC for short period collision) and a second 
one involving the third star (hereafter, LPC for long period col- 
lision). We can easily check whether any putative non-thermal 



800 



i 1 
i 1 
i 1 


1 1 1 1 1 
1 

1 

1 


' 1 
1 

1^ 

1 
1 




1 uu 

BOO 


i 1 

_ i i ^ 

- \ s ^ 


1 

i 
i 


1 
i 
i 
i 




600 


i i 
— i 1/^ 
i 


i 
i 
i 
i 


i 
i 
i 
i 




400 


i i 

yi i 

~ X i i 




i 




200 


y i ...--.--i-- 


i 


1 
i 


T 





5 10 15 20 25 

Wavelength (cm) 

Fig. 11. Dependence of the radius of the r = 1 radio photosphere 
as a function of wavelength for the primary (P), the secondary 
(S) and the tertiary (T) components. The vertical lines are lo- 
cated at 3.6, 6, 13, and 20 cm corresponding to the wavelengths 
investigated bv lBenaglia et al ] (l2006i) . We stress that the separa- 
tion between P and S is estimated at a sin ; = 28.9 whilst the 
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radio emission from the SPC could contribute to the measured 
non-thermal radio flux. The curves in Fig. [TT] are dominated by 
the primary contribution, which is easily explained by its high 
mass-loss rate. The radii at 3.6 and 6 cm are of the order of 
260 and 370/?© for the primary. The projected stellar separa- 
tion (ai sin / H- a2 sin /) is of about 28.9 which translates into a 
stellar separation of the order of 38 for an inclination of about 
49° (see Sect. l8.1l l. The radio photospheres at 3.6 and 6 cm over- 
whelm thus completely the close primary H- secondary system, 
and consequently the SPC is deeply buried in the opaque wind. 
This fact clearly rejects the possibility that the non-thermal ra- 
dio photons come from the interaction between the primary and 
secondary components. Alternatively, they could come from the 
LPC. The radii of the radio photosphere at 6 cm constitute there- 
fore minimum boundaries for the distance between the LPC and, 
respectively, the primary (solid curve in Fig. fTTT i and the ter- 
tiary (long-dashed curve). The sum of these two radii (i.e. about 
390 7?0) should be considered as a lower limit for the stellar 
separation between these two stars. That value is indeed much 
lower than the stellar separation estimated in Sect. 18.21 on the 
basis of the probable period of the wide system. A more strin- 
gent lower limit could be derived using values at 20 cm, since 
the free-free absorption increases as a function of wavelength. 
However, the non-thermal nature of the radio photons at that 
wavelength is not guaranteed by the spectral index reported by 
Benaglia et al. (2006). We note that the radio photospheric radii 
estimated here intimately depend on the assumed stellar wind 
parameters, which should be considered with caution. In par- 
ticular, the mass loss rate of the primary may be significantly 
underestimated here. However, our conclusion relies on conser- 
vative margins. 

9. Conclusion 

The analysis of high-resolution, high signal-to-noise data, col- 
lected over several epochs with FEROS, has allowed us to detect 
a third component physically bound to the inner system for the 
first time in the spectrum of HD 150136. In addition, we have 
revised the spectral classification of the three components to an 
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03V((f ))-3.5V((n) primary, an 05.5-6V((f)) secondary and 
an 06.5-7V((f)) tertiary. We also estimated their evolutionary 
masses as equal to 71, 36, and 27 Mq, respectively. 

We refined the orbital solution of the inner binary system 
by fitting the radial velocities of the different components of the 
SB3 system together The results have impro ved the accuracy of 
the m ajority of the parameters reported by Niem ela & Gamed 
(1200 5). Moreover, we found that the inner system is probably in 
synchronous co-rotation. The agreement between the evolution- 
ary and the Keplerian masses would suggest a probable inclina- 
tion of about 49°, implying dynamical masses of about 64, 40, 
and 35 Mq. 

From the systemic velocities measured on each epoch, we 
estimated that the orbit of the third star should have a period 
between 2950 and 5500 d and an eccentricity higher than 0.3. 
More data are necessary, however, to achieve an accurate orbital 
solution of the outer orbit. A possible alternative would be to 
observe HD 150136 by interferometry to constrain the missing 
orbital parameters. 

The determination of the individual stellar parameters of the 
three stars has yielded a luminosity and a gravity correspond- 
ing to main-sequence stars for the primary and the secondary 
components, but these parameters are trickier for the third star. 
Indeed, we observed a clear discrepancy between its luminosity 
(comparable to main-sequence stars) and its \ogg (typical of gi- 
ant stars). The origin of such a difference is thus unclear, but it 
seems that it could come from the disentangling procedure, no- 
tably because the spectral separation is not sufficient to sample 
the full width of the broadest lines, such as the Balmer lines, 
for the three components. Consequently, HD 150136 is probably 
composed of at least three 0-type stars located on the main se- 
quence and with inferred ages between 0-3 Myr. Finally, we 
showed that the non-thermal radio emission should come from 
the wind-wind interaction region between the third star and the 
inner short-period system, region that should be the site for par- 
ticule acceleration in this multiple system. 
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